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A SURVEY ON THE SPECIAL FUNCTION 



ANDREA OSSICINI 




Abstract. The purpose of the work is to furnish a complete study of a discrete and special 
function, discovered by the author and named with the Arabian letter ^ (Shinj^]. 

It includes three other papers, published in the international journal "Kragujevac 
Journal of Mathematics" . 

The methods, the techniques and the style of the demonstration inside such a work, are 
all related to Leonhard Euler, the very great Swiss mathematician. 



1 The letter ^ is the thirteenth letter of the Arabian alphabet. 
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Abstract. The purpose of the work is to furnish a study of a discrete and special function, 
discovered by the author and named with the Arabian letter ^J, (Shinjl, produced by a 

family of functions. A fundamental theorem is enunciated: it connects the rational values 
of this family with a natural number; to this aim two rational- value functions will be created, 
with the characteristic, in the field of real positive numbers, to be piecewise continuous. We 
state not only the existence of a separation element, but we prove that this element is 
just formed by one only integer constant function, the value of which is equal to 2. We 
point out a hypothetical, subtle connection among the special function j3 , the Eulerian 

function Gamma and the second-order Eulerian numbers. It is finally proved that ^ is 

completely monotonic: this characteristic is peculiar for the functions that have considerable 
applications in different fields of pure and applied Mathematics. 
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1. FROM THE CONTINUOUS TO THE DISCRETE 
Let's consider the three following transcendental functions, determined by the letter : 

. . 2fc+l 2fc+l •JfcJ-i 

us[i] = (i + ss) ;^[ 2 ] = ( 1 + 3^t) ;^[ 3 ] = ( 1 + 3^2 



and where the variable k will be restricted only in the real positive values, except some 
values, explained later. 

The analytic study of the three functions permits to identify their progress, for the 
presence, in each of them, of two asymptotes: a horizontal asymptote, got by the limit for 
k going to infinity, represented by the horizontal straight line, of heigho 

e 2/3 ps 1,947734041 
A vertical asymptote, specific for every function, got by the following limits: 



lim , £ [11 = oo ; lim , £ [21 = oo ; lim 



m + J [3] = 



OO 



Besides by the calculus of the first derivative of each function and the study of its sign, 
it's possible to verify that the ^£ functions are decreasing in their whole field of existence, 
more precisely, if we consider for the variable k the whole positive real axis, we must exclude 
at least for the second and third function, respectively the intervals ^0, I/3 and ^0, ^3 . 

In fact, if we derive the ^ functions we get respectively: 



_d_ 

dk 



2k+l 



j>[1] =J[1] •[2-log(l + i)-^ Tiy 
ij[2] =J[2] .[2.1og(l + 3^_) 
£ji[3] =J[3] ■ [2- log (l + g^j) 



fc(3fc-l) 
3-(2fc+l) 



-< V k G R + 

-< V k G R+ and k > V3 
-< V k G R+ and k> 2 / 3 



(3fc-l)-(3fc-2) 

In the Fig. 1 we have represented three functions where, among other things, it's evident 
that all of them have only one point of intersection with the horizontal straight line of height 
2. 



3 e represents the Euler's number. 
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horizontal asymptote 



1 4 C 



Fig.l 



That being stated, let's proceed in the passage from the continuous to the discrete, by 
considering for k, the only integer positive values; under these hypotheses it's possible to 
verify: 

jS [1] >- 2 for 1 < k < 8 and Jt [1] -< 2 for k > 9 
jSp] ^2for9</c<16 and ji [2] -< 2 for A; > 17 
jS[3] >- 2forl7</c<25 and j4[3] -< 2 for k > 26 

and therefore, in the discrete it's possible to define some interval^ Ig of integer values of 
the variable k, in order to characterize some limitations of the values that have the above 
stated functions; in fact it's possible to verify that for appropriate intervals If, £ = 2,3 the 
following are valid (see Fig. 2): 

ji [1] -< 2 -< j& [2] for 9 < k < 16 J 2 ; jS [2] -< 2 -< jS [3] for 17 < k < 25 => I 3 



4 £ shows the "ordinal number" of the interval h ■ 
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Fig. 2 



In conclusion it's possible to introduce a family of jS functions, by the definition of 
appropriate arcs, whose separation element is the horizontal straight line of height 2. 

The construction-algorithm of the above stated family is therefore describable in the 
following way: we begin from the first algebraic expression of the ji function, that is 

/ 1 \ 2k+l 

( 1 + 3^ J , starting to calculate by growing values of the integer positive variable k, 
the corresponding rational values of ; for the first 8 integer values of k, the function 

has rational values greater than 2 and it's therefore possible to associate to such values a 
bounded arc of the same function, represented in the discrete field by a sequence of rational 
numbers, each of them greater than 2. 

After that we decrease of a unity the value of the denominator of the fraction inside 
the ji function, consequently we'll get the algebraic expression of a new function, that is 

( 1 + 3F=T ) ' * na * we nave previously identified with ^ [2] and that we can define as 
the following of ^ [1] . 
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We'll repeat, for it too, the same procedure and therefore we'll calculate by growing 
values of the integer variable k, but greater than the previous 8, the corresponding rational 
values of [2] ; also in this case for 8 integer values of k , the function has rational values 
greater than 2 and therefore it's possible to associate to such values a bounded arc of the 
same function, represented in the discrete field by a sequence of rational numbers, everyone 
greater than 2. 

Now, for the same 8 integer values of k, utilized for ^ [2] , it's besides possible to verify 

that the ^ [1] function has, on the contrary, rational values smaller than 2 and it's therefore 
possible to associate to such values a bounded arc of the same function, represented in the 
discrete field by a sequence of rational numbers, all smaller than 2. 

If we repeat the procedure and then we decrease, as usual, the value of the denominator 
of the fraction inside the ^ function, we can build the ^J, [3] function, that so results the 

following of ji[2] . 

In this case, differently from the first two, exactly for 9 integer values of k, greater than 
the previous 8, the function has rational values greater than 2 and therefore it's possible 
to associate to such values a bounded arc of the same function, represented in the discrete 
field by a sequence of rational numbers, all of them greater than 2. 

Here, for the same 9 integer values of k, it's possible to verify that the j3 [2] function 
has, on the contrary, rational values smaller than 2, too, and it's so possible to associate 
to such values a bounded arc of the same function, represented, in the discrete field, by a 
sequence of rational numbers, each of them smaller than 2. 

Consequently about what described, if we consider the second interval of 8 integer values 
of k, that is k =9,10,11,12,13,14,15,16, we can build two arcs, represented by two sequences 
of rational numbers: the first sequence formed by numbers greater than 2, because belonging 
to the arc of the <ji[2] function, the second sequence formed by numbers smaller than 2, 

because belonging to the arc of the <ji[l] function (see Fig. 2). 

What shown is repeatable and it's possible to experiment, while the integer variable k 
grows, the determination of two appropriate arcs, belonging to two following ^ functions. 

To simplify the use and the control of the described algorithm a vector function has 
been defined by the software product DERIVI^ Version 6 for WINDOWS, which allows 
the display of 11 consecutive values of a generic ^ function, by two only parameters: an 
integer value of k and a further integer value, corresponding to the value of the interval 
(decreased of a unity) that we wish to study. 

In APPENDIX it's given, besides the macro function, which identifies the vector func- 
tion, the result of a display got by its use. 



5 DERIVE is a powerful instrument of CAS (Computer Algebra System ), spread by Texas 
Instruments. 
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2. THE FUNDAMENTAL THEOREM OF THE , i FUNCTION 



Generalizing what described in the first paragraph and using the same method to build a 
couple of functions, identifiable by a precise interval 1%, it's possible therefore to enunciate 
the following Fundamental Theorem: 

Let's consider k,£ natural numbers different from zero, an interval of integer values 
of the variable k, $1(1 t) an auxiliary integer function, then for each interval It it's always 
possible the construction of a couple of functions with rational values, exclusively depending 
on k and that we denote with i Ji, so as the following boundary is always valia®. 

j[k,n (it - 1)] -< 2 -< j>[k,n(i e )] with o (h-i) = n (i £ ) - 1 ; vi t , mgn 

and for the generic ^ function it 's valid: 

/ i \ 2fc+l 

^"■ n ^ = ( 1 + 3^TT(7ri) « 

The auxiliary integer function $l(Ie), that really represents a growing "step function", 
is defined, for the intervals of 8 or 9 following values of k, in the following way: 



• 


n 


(h) 


= for k=l,. . . 


,8 






• 


n 


(h) 


= 1 for k=9,. . . 


,16 ; $1 (J 3 )=2 for k=l7,. . 


. ,25 ; fl (h)=3 for k=26,. . . 


,34 


• 


n 


(h) 


= 4 for fc=35,. . 


. ,43 ; Q (I fl ) =5 for k=A4, 


. . . ,51 ; n (J 7 )=6 for fe=52,. 


. . ,60 


• 


n 


ih) 


= 7 for k=61,. . 


. ,69; n (I 9 ) =8 for k=70,. 


..,78 ; Q (Ji )=9 for k=79,. 


..,86 


• 


n 


(hi. 


) = 10 for A;=87, 


. ..,95 ; $1 (7 12 )=11 for k= 


96,. . . ,104. ; etc. 





The sequence of growing values, defined for the auxiliary [^function $l(Ie), is therefore 
obtained to allow the individuation of two bounded arcs, belonging to two following ^Jb 
functions, one above the line of height 2 and the other below it. 

The set of the 11 specific intervals, by which the function O(I^) is always positive and 
growing, represents a dominant characteristic of the family of ^ functions. 

For a better precision we write the numerical series that identifies the extent of such 
intervals, in terms of consecutive values of the variable k: 



8, 9, 9, 9, 8, 9, 9, 9, 8, 9, 9 for a sum of 96 values of the variable k. 



6 The boundary can include the sign "=" if the integer variable k goes towards the infinity. 

7 The function O(J^) is a depending function on the integer variable k. 
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As we can see, the first interval of 8 values has been neglected: in it the function f2(/i) 
has value zero, but such interval has a peculiarity, as exactly at interval 31 the above 
mentioned series is interrupted and an interval of 8 values takes placed of one of 9, giving 
origin again to the same sequence: the following 11 values of such function are determined 
by the extent of intervals, typical of the series. 

Now this event is regularly repeated and precisely every 40 and 51 intervals: for instance 
the first substitutions happen at the intervals: 31, 71, 122, 162, 213, 253, 293, 344, 384, 435, 
475, 526,566,617,657,697,748,788,839,879,930,970,1021,1061,1112,1152 and, all of them, in 
advantage of the 8 value intervals. 

The first real effect of this phenomenon on a procedure that allows to calculate exactly 
the last value of k, present in a determined interval, foreseeing a constant and complete 
repetition of the numerical series of 11 intervals, is exactly quantifiable at the interval 122, 
where it is practically possible to verify that the last value of k inside it, results lower of 
a unity: this means that we could obtain the same result if at the interval 122, excluding 
the initial one, there were 120 integer intervals and one reduced of a unity, because of one 
only interruption. 

Growing k, the substitutions immediately determinate a further but stable effect, for 
instance, at the interval 617 it's possible to verify that the last value of k, present in the 
interval is of two unities inferior than the computable one and, going further it's at the 
interval 1112, that it is possible to verify that the last value of k, present in the interval, is 
of three unities inferior than the computable one, and nothing short, for many following 
intervals it's possible to observe that such reduction of a unity is noticeable exactly every 
495 intervals for indeed 10 times (first case) and every 484 intervals only once (second 
case). 

In the first case we can get the same results by hypothesizing the constant presence of 
494 integral intervals and one, reduced of a unity and this because 11 substitutions produce 
inside 495 intervals 38 complete series of 96 values of k and 11 series, reduced to only 61 
values of k; under these conditions it's immediate to verify that we obtain 4319 values of 
k, that are inferior of one only unity in relation with the possible values (4320) inside 495 
intervals, in case there were exclusively 45 numeric complete series. 

In the second case we can get the same results by assuming the constant presence of 
483 integral intervals and one, reduced of a unity and this because 11 substitutions produce 
inside 484 intervals 37 complete series of 96 values of k and 11 series, reduced to only 61 
values of k; under these conditions it's immediate to verify that we obtain 4223 values of 
k, that are inferior of one only unity in relation with the possible values (4224) inside 484 
intervals, in case there were exclusively 44 numeric complete series. 



8 For precision it is the last 9, belonging to the second group of three consecutive 9 , that is of 
the eighth term of the series. 
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The peculiarity of the first interval of 8 values, for which the function O(I^) is worthless, 
consists therefore in the fact that, as it was determined before the birth of the sequence, 
characteristic of the series, it is possible to think that at the origin, the number of integer 
intervals corresponds to 1. 

We specify that such characteristics, even if relevant, are minor, both in relation to the 
continuous and repetitive presence of the above described series and in relation with the 
largeness of the single intervals, which never descend under 8 values^ for very great values 
of k. 

The ratio of the number of integers inside the numeric series with the 11 intervals 
corresponds to 8, 72 and this value diminishes in a little meaningful way while k tends to 
infinity, if we consider the ratio of the value of a very great k with its own belonging interval 
If. a sufficiently precise value is obtainable with the following expression 



96 



1 



TT 



where 



we deduce by the denominator of the second fraction the importance of the number 494. 

Successively we give (Fig. 3) the graphs, related both to the family of ^Ji functions, or 
better, to the set of the arcs belonging to them, and to the auxiliary function and suc- 
cessively in particular the development of the graph of a couple of functions, characteristic 
of a precise interval Ii (Fig. 4) and of the pointers than put into evidence the behaviour, 
growing k. 

This last behaviour results particularly evident, by examining the various displays, pro- 
duced by DERIVE, in relation with the first 11 intervals, typical of the standard series, 
while for the greater values of k it's necessary to outdistance in an appropriate way the 
intervals on which to do a comparison to have a further confirmation of such behaviour. 

Now if k tends to infinity it's possible to compute the limit towards which the generic 
j2 function, that in this case represents the fusion, to infinity, of two arcs of following ^ 
functions, that refer themselves to a hypothetical and extreme interval In. 

II- 

115 

■I 

II- 
■■' 

■:. - 
4- 



3- 

: - 

i - 



"0 s |*,na-i)j 



"Tr- 
ie 



~n n — 

It 1 7 is-n 



~n n — 

u is n u 



~n — 

■51 !J 



Tl 



"Tl 

mi Til 



Tl n — 



TF i — 

IM 



Fig. 3 



9 This means that every interval It will never be empty. 
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Keeping in mind what noticed in the first and second case, previously described we can, 
first of all, calculate the two following values, by "excess" and "defect" of the ratio of a 
value of k, with the value of the integer function 0,(lt), determined as a function of the It 
interval, containing the same value of k: 



I st case: 



A: 



4319 



n(I £ ) 495 



8,725 



2 nd case:- 



k 



4223 



8, 7252066115702479333884297 



n(I e ) ' 484 

These values, by considering the frequency of the two cases ( 10 times the first case and 
once the second case), permit to the approximate the ratio with: 



k 



10 * 8, 725 + 8, 7252066115702479333884297 



11 



= 8, 7252483512814091326488 



96 



1 



96 
IT 



1 

494 



11 493,9793814432989690721649 
In this way, according to the analysis of the progress of the various and following 

functions, we can furnish only an esteem of the ratio n . 

If we want to determinate the effective value of such ratio, it's necessary to compare the 
values of the two quantities for sufficiently great values of k. 

Keeping in mind the algorithm described at the end of paragraph 1, which puts in 
evidence the continuous oscillation of the rational values of the following ^ functions, in 

proximity of the integer number 2, we will compare the values of the ratio n ^ as k grows 
with the value: (3 - 2/log2) _1 . 




Fig. 4 
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More precisely, considering the values raised to the 18 power (ki = 10 18 , that is one 
quintillion) , raised to the 33 power (k 2 = 10 33 , that is one decillion) and indeed raised 
to the 63 power (k 3 = 10 63 , that is one vigintillion) it's possible to observe that the real 
value of the ratio ^ (k) = n jj ^ tends to the previously shown value; in fact in the three 
various cases it's possible to determinate in order what follows: 

ttfa) = TjTTT - (3 - 2/log 2)- 1 -< 5 • 10~ 19 ; tt(A&) - (3 - 2/log ~< 5 • 10" 34 ; 

=► *(Jfea) - (3 - 2/log 2H 1 -< 5 • lO" 64 
Definitively, going to the limit for k, which tends to infinity we can expect that: 

lim (-pr^-) = (3- 2/log 2)- 1 



This result has obviously an immediate consequence in the calculation of the following 
limit: 

i \ 2fc+l 



lim J[k,n (It)] = lim 1 + 



k^oo \ 3k — Q (Ii) 



■2j, i l Hnu 2^±L. 



Zk-\-l lili - L fc — too / 9 \ 

lim e 3fc-n(/ £ ) =e 3k - k { 3 -T^2) =e lo § 2 



In conclusion, it happened that the natural number 2 (by the integer constant function 
Y=2) can represent the separation element of two arcs, belonging to two following ^ 

functions, which therefore can be identified with two contiguous classes^!, which tend to 
approach indefinitely. 

Even if we observe that the validity of the fundamental theorem is actually included 
in the construction algorithnj^l of the family of functions, keeping obviously in mind the 

characteristic of monotonicity (see previous paragraph and the following) of the generic 
function, by a more appropriate notation, due to Iverson, further we will give the elements 
to get a strict proof of the theorem in the modern sense of the term. 

Let's extend the dependence of the integer step function £1 (//) to the real field and let's 
use the following definition: 

I \ 2x+l 



(x) = min {k G N: Sfc + i (x) > 2 } ; x G R and where (x) = 1 + 



3s- Jfe + 1, 



10 The contiguous classes are meant represented by two groups of rational non- integer numbers, 
greater and smaller than 2, separated by the rational integer number 2, which obviously doesn't 
belong to any of the two classes. 

11 Let's observe that the algorithm determinates for the generic function an inferior limit: it's 
established an inferior extreme which corresponds to the integer number 2. 
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By simple algebraic passages we have that: 



n (x) 



3x 



22^+1 - i 



(2) 



where [~x] means the smallest integer, greater than x or equal to it. 

So, from such a particular definition of f2 (x) the fundamental theorem immediately 
derives. 



In fact the ^Jb function possesses the following explicit formula: 

( 



S(x) 



1 + 



1 



22+1 



3x — (x) 
and passing to the sequence we have: 

S(n) = 



1 + 



V 



3x — 


3x ji 







2x+l 



1 + 



V 



2 2n + l _1 



2n+l 



where [xj shows the greatest integer smaller or equal to x. 

From now on we will use the locution of " special function ^ " for the discrete function 

^ (I^)] , defined by the formula (p}. 
Extending the field of definition of the variable k to the real positive numbers, it's 
possible to notice that such function, being represented by the union of continuous arcs 
(all above the straight line of height 2 ) is actually assimilable to a piecewise continuous 
function. 



3. AN APPROACH OF THE SPECIAL FUNCTION Js WITH THE EULERIAN^ 

GAMMA FUNCTION T 



One of the fundamental characteristics of the Eulerian gamma function is a certain 
condition of monotonicity, which is the fact that the function, which can be intended as the 
most spontaneous extension of the factorial n \ = 1 ■ 2 ■ 3 ■ ■ ■ n, out the field of the natural 
numbers, is logarithmically convex. 



Chapter II, [5]. 
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If we consider the expression (prj and we calculate its first derivative, we have: 

d 



dk 



j[k,n (h) 



ji [k, ft (It 



V k G N and & > 



2 log 1 + 



1 



3- (2k + 1) 



(3k - ft (It)) J (3k - ft (I e )) (3k - ft (I e ) + 1) 



(3) 



n (ia 



it is always negative, therefore the function results monotonic and 
decreasing, and besides it's interesting to observe that, growing the k value, its absolute 
value diminishes. 

That being stated we can also verify that the special function has the same charac- 
teristic of monotonicity of the Eulerian function T: it's sufficient in fact to verify that the 
second derivative of the logarithm of the same j3 function is positive. 

Such result is immediate, in fact for the ([3]), being: 

^\o g j[k,n(i t )} =A ( ±j[k,ci &)] / j[k,n m > 



we have that: 



d_ 

dk 



21og 1 + 



1 



(3k - n(i e )) 

12 • (3k ■ n (It) - 



3- (2k + 1) 



(3jfe - n(i t )) ■ (3k - n(i e ) + 1) 
n (h) 2 ) + 6 • (6k - n (h)) + 9 



(3k - ft (l e )) 2 -(3k-ft (I £ ) + l) 2 
Since, for definition we have k >~ ft (It) ,V k E N, we have in conclusion: 

^o g j[k,n(i £ )] y o 

This demonstrates that also the special function ^Ji is logarithmically convex. 

The special function ^ is actually an exponential general function; that being stated, 
considering its base, defined by the analysis of its behaviour to infinity, we consider the 
following two particular expressions: 



n(i+ 

k=l 



log 2 
2k 



J and 2_j 

' fc=i 



log 2 
2k + log 2 



In the previous paragraph we have practically put into evidence that for great values of 
k the rational non-integer term of the base of the special function tends (~) to: 



1 log 2 

3k - ft (I e ) 2k~ 



(4) 
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consequently, as there is the following relation O ( Ii — 1) = O (Jg) — 1, we have: 

1 log 2 



(5) 



3fc-fi(/f-l) 2/c + log2 

The value (JH), if fixed for the base of the generic ^3 function, is very well fit to "inter- 
polate" in the continuous field the union of the various arcs of the family of jS functions, 
got in the discrete field, above the straight line of height 2. 

By this last value it's possible to calculate and verif\f^l the following sizeable expression: 



l^V 2k ) r(n + i).r(i + '- 2X 



The value ©, on the contrary, if fixed for the base of the generic ^ji function, is very 
well fit to "interpolate" in the continuous field the union of the various arcs of the family 
of ^ functions, got in the discrete field, below the straight line of height 2. 

By it, it's possible to define the following partial sum: 

" log 2 log 2 ™ 1 m 

^2* + log 2 2 f^k + ^f U 

From here, keeping in mind the formulae of recurrence of the logarithmic derivative of 
the gamma, function named digamma : 

ip(x + 1) = ib{x) H — ; ibix + n) = ib(x) H 1 1 h H with n > 1 

x x x+1 x+2 x+n—1 

we have, giving x the value 1 + ^f-^, the following identity: 



log 2 _ log 2 
^ 2fc + log 2 2 

Now, going to the limit for n — > oo, for the known properties of the digamma function, 
we have that the following series is divergent, that is: 

^ log 2 

That being stated, the importance of the result ([6]) must be evaluated above all according 
to the following identity: 

1 + log2\ log 2 



2k J V 2/c + log2 



13 



r (x + n) = x ■ (x + 1) • (x + 2) • • • (x + n - 1) • T (x) with n G N and x !>- 
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From this last identity we can deduce: 



k=l 



log 2 \ r (n + 1) • r (1 + 



TT i _ ^ " = —— V L_f J (g) 



and considering the infinite product: 



k=l 



log 2 
2k + log 2 



diverges to zero. 

Such conclusion can be immediately verified also by applying the result dH|), modified in 
the sign and considering the following: 

Theorem: Supposed that —1 -< a n < and the series >~2 a n is divergent, then the 
infinite product FJ (1 + a n ) diverges to zero. 

Proof 

Assumed b n = — a n =>■ < 6„ -< 1 , as for < x -< 1 results 1 — x < e' x , we can 

write: 

n 

-< P n = J] (1 + or) < e-( bl+b2+ +K) (10) 
r=l 

then if the series X] a n is not convergent and necessarily diverges to — oo, from (|10p we have 
that P n —* that is the infinite product diverges to zero. 

Besides, keeping in mind that, supposed a and h non negative, the following relation is 
valid: 

^ r • (r + q + b) = T (1 + a) • T (1 + b) 
1-1; ( r + o)-(r + 6) ~ T(l + a + 6) 

verifiable, applying the known formula by Euler: 

r fx) = lim 



rwoo a; • (a; + 1) • (z + 2) • • • • (s + n) 
and considering the limit for n — > oo of the development of the following finished product: 



n 



r ■ (r + a + b) _ (1 + a + 6) • (2 + a + 6) • • • (1 + a + b + n) 



L (r + a) • (r + b) (n 1+a+b ) • n! 

(n 1+a ) • n! (n 1+fe ) • n! n + 1 



(l + a)-(2 + a)---(l + a + n) (1 + 6) • (2 + b) ■ ■ • (1 + b + n) n 



Chapter II, pag. 33, [6]. 
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we can state that: 



r(l + n) • T (l + i2|2) °° r-(r + n + 



n 



log 2 
2 



r (l + n + ^ (r + n) ■ (r + ^ 

By exploiting the properties of the gamma function, or resorting to the definition of the 
famous hypergeometric series, it's also possible to verify the following: 

r(i + n)-r(i + ^) _^ (-ir-(^)-(^-i)---(^-r) 
r (i + n+ io|2^ A. (1 + n + r ). r! 

In fact, adopting the classic symbolism for the hypergeometric functiorl^l : 
„, , s a-b a ■ (a + 1) • b- (b+ 1) , 

a ■ (a + 1) • (q + 2) • b ■ (b + 1) • (b + 2) 3 
+ c-(c + l)-(c+2)-3! X + "" 

we can exploit the possibility to express it in terms of gamma functions, considering the 
hypergeometric integral, that is: 

l 

F (a ' 6; C; X) = T(b) T T(c-b) ' / tb ~ K (1 " t)C ~ b ' 1 ' (1 " 1 ' ^ ^ (5?C ^ m ^ 0) 

From this, considering the limit for x — > 1~ (Abel's theorem) and exploiting the proper- 
ties of the Eulerian Beta function, we obtain the important relation of the hypergeometric 
Gauss's series: 

F(a,b;c; 1) = ^ • r . (c ~ Q ~^ (c# 0,-1, -2,..., K(c-a-6>0) 

r (c — a) • r (c — oj 

With the positions a = — -^7p, b = n and c = n + 1, as all the required limitations are 
satisfied for the parameters, we easily reach the (jll|) . 

In conclusion, starting from the characteristic base of the generic jS function, defined by 
the analysis of its behaviour to infinity, we have also stated the following sizeable relation 
between an infinite product and a numeric series: 

( log 2 
F( — 1_ n; n + l; i 

oo r-(r + n + ^) _ ~ (-1)" ■ • - l) • • • - r) 

11 / \ / . 1dp-2\ ^ 



= 1 (r + n) (r + ^) £5 (1 + n + r) 



In the real field such result is graphically represented by Fig. 5, shown below, where it 
is evident that the x axis represents a horizontal asymptote for x — ► +oo. 



15 Chapter III, [5]. 
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Besides, in the field of the non- negative real numbers, both the infinite product and the 
series are absolutely and uniformly convergent. 
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Fig. 5 



4. THE SPECIAL FUNCTION J AND THE SECOND-ORDER EULERIAN 

NUMBERS^ 

In the second paragraph we have shown that one of the most important characteristics 
in the construction of the family of ^ functions is related with the repetitive presence of 
11 specific intervals Ig of the integer variable k: we have indeed seen that the extension of 
such 11 intervals is always characterized by the following numeric series: 

8, 9, 9, 9, 8, 9, 9, 9, 8, 9, 9 for a sum of 96 values of the variable k. 

But it also happened that such series undergoes inside a precise number of integer 
intervals some interruptions, as 8 value intervals substitute some 9. 

More precisely this event is regularly repeated every 40 and 51 intervals and it's possible 
to observe that such intervals are aggregable so that to give origin to two groups: the first 
one formed by 484 intervals, produced by 7 interruptions at a distance of 40 intervals 
and 4 interruptions at a distance of 51 intervals, the second one formed by 495 intervals, 
produced by 6 interruptions at a distance of 40 intervals and 5 interruptions at a distance 
of 51 intervals. 



pag. 247-251, [3]. 
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For each group of 484 intervals we besides note always 10 groups of 495 intervals: this 
corresponds on average to a number of 494 intervals for all the groups. 

In conclusion the characteristics of such phenomenon, keeping also in mind what ex- 
plained in the initial part of the second paragraph, allow correctly some estimates for the 
integer values of k, referred to an interval Ig, if we suppose the effective presence, even if 
virtual, of the following sequence of integer intervals: 

1,120,494,494,. . . ,494, ^4943,494,.., 494,. . . ,494, ..,494,. . . , «494 (ad infinitum). 

This interpretation is besides confirmed by the check carried out on the limit of the 
ratio % fay 

Now, as these numbers are closely connected to the nature of the complete and periodic 
sequence of the numerical series shown above, we have supposed that the numbers 1 , 120 
and 494 can belong to a category of special numbers; particularly a research pointed to 
verify this hypothesis, has implied the following curious discovery even if it is always a 
conjecture on an almost light and remote connection: the numbers 1, 120 and 494 belong 
to the family of the so called "second-order Eulerian numbers" , which result important for 
the tight connection that they have with Stirling's numbers. 

They satisfy similar recurrence to the characteristic one of the "ordinary Eulerian num- 
bers" , which are useful, above all because they give a connection among ordinary powers 
and consecutive binomial coefficients. 

To be clearer we show below the recurrence, which characterizes the second-order Eule- 
rian numbers, showing how it's possible, by it, to produce the three special numbers, typical 
of the special function ji. 



In fact we have: 




(k + 1 




+ (2n - 1 





1 Vn/ 




for n 7^ 



17 «494: with such notation we want to put into evidence that for great values of k the estimable 
number of integer intervals is, in some very near case, to such value, in fact it's as if every 495.000 
intervals, two others of them, on average, undergo a reduction of one unity (from 9 to 8), but on the 
whole this further phenomenon is absolutely neglectable. 
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6 \ \ / / 5 \ \ , / / 5 



1 



(2 > +(10 > o = 114 



I -» ! » o H 4 ° 



8 \\ // 7 \\ „ JN // 7 



(2) i + (14) o = 494 



At the end, to be complete, we give besides the general formula of the second-order 
Eulerian numbers, which puts into evidence the connection with the binomial coefficients 
and with Stirling's numbers, a representation of the second order Eulerian triangle (Fig. 6): 



((:)>- 1 ( y 1 ) { n :+TV k } ■ ^ * » ° <-> 

where for the binomial coefficients it is valid: 



k 1 - k ! (n-fc) ! 
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and for the 
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Fig. 6: Second Order Eulerian Triangle 

If we consider, for fundamental of the special function the number / / ^ \ \ from 
(fT2l) we have: 

«?M")("h(")(i) 

but keeping in mind that for Stirling's numbers are valuable the following identities for 
n y 0: 

n 



we finally have: 



2 [ = 2 - 1 and 1 =1 



1 • ( 2 9 - 1 ) - 17 • 1 = 511 - 17 = 494 



5. FROM DISCRETE TO COMPLEX FIELD 



In the previous paragraphs we have discussed on arguments by infinitesimal, asymptotic, 
numerical and combinatorial analysis to characterize the special function ^ . 

It is known that in these fields the completely monotonic functions play a fundamental 
role. 

We recall that a function / : / — > 5ft is said to be completely monotonic (cm.) on a real 
interval I, if / has derivatives of all orders on / which alternate successively in sign, that 
is: 

. /W ( x ) > Vx e 7 and V n > with n = 0, 1, 2, 3, ... 
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In the recent past, various authors^ showed that numerous functions, which are defined 
in terms of gamma, polygamma and other special functions, as the hyper geometrical ones, 
are completely monotonic and used this fact to derive many interesting new inequalities. 

We shall confine ourselves only to prove, in very simple way, that the special function ^Jb 
too, even if it is for definition a piecewise continuous function, in the real field, it possesses 
the same property (cm.). 

Lemma. If f (x) and g (x) are cm., then a ■ f (x) + b ■ g (x) , where a and b are 
non-negative constants and f (x) ■ g (x) are also cm.. 

The proof of the first thesis is obvious, the second one is then easily seen from the 
Leibniz formula : 



dx n w v " y 3 V ~ /J ^ \ k 



f{x) -g (x)] = £ ? f (k) (x) ■ g^ (x) 



Theorem 1. The special function j5 [k,Q (Ie)] is completely monotonic in each In . 

Proof: The base function / (k) = 1 + 3fc _^ ^ is cm. in each Ig ; in fact the n-th 
derivative of this function is: 

^ /(fc > ^» - (-i)" . /<») > o 

Thus by the Lemma, in case k is an integer, it's obvious that also the special function 
^ is cm. in each Ip ; therefore it remains to prove Theorem 1 for the case when k is a real 

and positive numbeJ^l. 

To obtain this it is necessary to the use the following obvious Theorem 2, which is a 
consequence of the Lemma, for composed functions. 

oo 

Theorem 2: Let y = f (x) cm. and let the power series ip{y) = a jU 3 converge for 

j=0 

all y in the range of the function y = f (x) . // aj > for all j = 0, 1, 2, 3, ... then ip [f (x)} 



is cm.. 



Corollary: // / (x) is cm., then e^ x ^ is cm.. 
In particular, as the special function ^J, is equal: 

j[k,nm = e (2fe+1) - los ( 1+ ^w) 

it is cm. in each If . 



18 pag. 445-460, [1]. 

19 We must also remember the ([2]) without giving up the characterization, determined for each 
interval Ie . 
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An interesting exposition of the main results on completely monotonic functions is given 
in Widder'a^l work. 

That being stated, with some limitations due to the nature of the special function 
in the real field, we can describe a method for estimating the same function in the complex 
field with an important improper integral. 

The ^ function results to be a piecewise continuous function because of the presence of 
the step function f2 (It), that is discontinuous, and actually its complete monotonicity has 
been proved for all the closed intervals Ii . 

Passing to the interval I [0, oo ) the continuity is not guaranteed and therefore the ap- 
plication of the following Hausdorff-Bernstein-Widder'a^] Theorem, must be done carefully, 
or it can be limited to characterize the behaviour of the special function ^ at the origin 
and infinity. 

Theorem 3: A necessary and sufficient condition for the function f (s) in order to be 
completely monotonic in the interval I [0, oo ) is that: 

oo 

f(s) =L s [F(t)}= Je- st dF(t) (13) 
o 

where F (t) is non- decreasing and the integral converges in the interval I [0, oo ) . 



The (I13p represents the transformation of Laplace-Stieltjes of a locally and absolutely 
continuous function, with real values, in the interval I [0, oo ) . 

Now, keeping in mind the relation existing in such case between Laplace-Stieltjes trans- 
form L, and „,*n,,n@ Laplace transform L < observing that we can suppose F (0) = 0): 

oo 

L s [F (*)] =s-j e~ st F (t)-dt = s-L[F (t)} (14) 
o 

we can determine the expression of the function F (t) by the inverse Laplace transform. 

From (|13|) and (|14p we get (s) = = L [F (t)] and successively we need to face the 
calculus of integral of the tvpd^l: 

I 

XQ—i-OG 



1 X0~\-l-OG 

F(t) =L- 1 [$(s)} =v.p.-± l - J e ts $(s)-ds (14.1) 



with t y- , (s) holomorphic function in the half-plane 5Re (s) >- and xq arbitrary real 
positive number. 



20 [?] 

21 Chapter IV, pag. 160-161, [7]. 

22 Chapter IV, [2]. 

23 pag. 168, [2]. 
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In general the calculus of Bromwich's integral (14.1), that is defined as a Cauchy 
principal value, we can only do it numerically, applying quadrature formulae, but in our case 
[ f (s) = jJ (s) ], being the (s) function piecewise analytic, we could prove by the direct 
calculus of the considered integral (see the following paragraph), that F (t) = 2-0 [t] and 

[t] represents the unit step function or Heaviside's function ( [ t < ] = ; [ t > ] = 

1 ) and in ([13]) we'll have dF (t) = 2 • d@ [t] = 2 • 5 (t) dt, with 5 (t) that is the distribution 
of DiracS 

6. THE BEHAVIOUR OF THE SPECIAL FUNCTION J AT THE ORIGIN AND 

INFINITES 

In the previous paragraph we have stated the following approximation "~" , in terms of 
Laplace transform: 

oo 

*> 00 = J («) / B ~ L [F (t)] = L [2 • (t)\ = J e~ st 2 ■ (t) dt = I 

o 

The approximation is essentially origined by neglecting the point of discontinuities of 
the first kind of the special function j3 , between an interval Ii and the following I^i as 
far as the interval / [0, oo ) . 

Said that to calculate the Bromwich's integral (14.1), we consider the path of integration 
rightly deformed, as we can see in the Fig. 7. 
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Fig. 7: The path of integration 



pag. 30-35, [2]. 
pag. 491-492, [4]. 
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<3? (s) is a multiple-valued analytic function; in the interest of obtaining a single value 
mapping we consider the principal branch of the power in (s) , that allows us to treat a 
holomorphic branch of the same $ (s) function. 

In its field of single- values, the function $ (s) possesses a simple pole at s = and two 
branch points at s = — | and s = — | ; let's observe, in fact, that in the interval [—1, 0] the 

special function J is: J (s) = (l + = e {2s+iyi ° s ( 1+ ^ 

A cut, joining the two branch points, would prevent s to circulate around them, and 
the special function ^ (s) can be treated as a piecewise holomorphic function. 

In conclusion, defined with C the boundary (shown in Fig. 7), with Cr the circular 
arc of radius R, with B the vertical line with 3ft (s) = xq, with L the boundary of the cut 
branch, formed by D £ and E £ , that is the semi-circles of radius e capping the ends of the 
branch cut and by L + and L~ , lines above and below, we'll have: 

xo+i-oo 

/ ^ W .* = JL{ /...-/...-/...-/...-/...-/...} (is) 

xo-i-oo C Cr L+ L- D e £ £ 

The first integral in the second member of (|15p . by the residue theorem, is 2 • O (i), as 
the first order pole, at the origin, gives: lim^o ^ (s) = 2. 

The second integral with s = R ■ e l ® is: 

tt/2 3?r/2 3tt/2+<5 

J e ts $(s) ■ ds = J ...d$+ J ...d#+ J ...d-d (16) 

C R 7r/2-<5 tt/2 3tt/2 

and therefore it vanishes: in fact the first and third integral in the second member of (|16j) 
vanish as R — > oo by the maximum modulus bound and the second integral vanished by 
Jordan's Lemma. 

The third and fourth integral in the second member of (|15p cancel each other along the 
paths L + and L~: in fact their values, calculated along their opposite paths, eliminate each 
other. 

In the end, the last two integrals of (|15p. by the maximum modulus bound, vanish as 
e -» 0. 

Now, from the asymptotic behaviour of the function F (t), so determined, we can there- 
fore deduce asymptotic properties of the correspondent Laplace transform, that is to use 
the following Abelian theorems (initial and final value theorem): 

Theorem: Let F be a transformable function and let's suppose that the ]hn.t—>ooF{t) 
exists, then the lim^o [ s ■ $ (s) ] exists, too (let's suppose, for convenience s £ di) and 
is: 

lim [ s ■ $ (s) ] = lim F (t) (17) 

s^O t— >oo 

// the lim t _ +0 + F (t) exists, then the lim^oo s ■ $ (s) exists, too and is (with s £ 3?) : 

lim [*•*(*) 1 = lim F (t) (18) 
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In our case we have seen that F (t) = 2 • (t) , with (i), that is Heaviside's function, 
and therefore we have that: 



and also: 



lim F (t) = lim 2 • G (t) = 2 

t— >oo t^oo 



lim F (t) = lim 2 • 9 (t) = 2 



The existence and the calculus of such limits, in (jlT[) and (I18p . give the following results 
(between them a further confirm of the fundamental theorem): 

lim . i (s) = 2 ; lim , i (s) =2 
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APPENDIX 

DERIVE Version 6 : the vector function : S(k,fi) := lim t ^. fc (1 + gj^J 
#1: V(k, n ) := VECTOR([t, o, S(t,Q )], t, k, k + 10) 



a display result: 



V( k, 


n 


) VECTOR 


09 


i 


2.0484148121729077984789748528464903812082646338028 


10 


i 


2.0379259208387064562838079920238964441117176933744 


11 


i 


2.0294161672236677191636908626945714916029532064944 


12 


i 


2.0223737073469397533461445484297949184415534016560 


13 


i 


2.0164491799135882361365114303236301480590762568587 


14 


i 


2.0113959747189663594458436566806886371655211482846 


15 


i 


2.0070350457364044054984268130457357906298862812402 


16 


i 


2.0032332566108411453651981972386971733090123101528 


17 


i 


1.9998895526624551656968593976078763119133058370586 


18 


i 


1.9969258468076576081148471529242828805292201418873 


19 


i 


1.9942808454379732420411582337304540085201659981375 



26 



Kragujevac J. Math. 29 (2006) 26-35. 



THE SPECIAL FUNCTION J , II. 



Andrea Ossicini 



Via delle Azzorre 352-D2, 00121 Roma, Italy. 
(e-mail: a.ossicini@finsiel.it) 

(submitted September 9, 2005) 

Abstract. We describe a method for estimating the special function ji , in the complex 
cut plane A = C\ (— oo,0], with a Stieltjes transform, which implies that the function 

is logarithmically completely monotonic. To be complete, we find a nearly exact integral 
representation. At the end, we also establish that l/ ( ji(x) is a complete Bernstein function 

and we give the representation formula which is analogous to the Levy-Khinchin formula. 
2000 Mathematics Subject Classification: primary 33B15; secondary 26A48; 33E20. 

Keywords: special functions; completely monotonic functions; integral transforms; Bern- 
stein functions. 

1. INTRODUCTION 

In [8] the author introduces a new special function, named with the Arabian lettei0 ji, 
and proves that this is logarithmically convex and completely monotonic for all the closed 
real intervals It with £ = 1,2, 3, ... . 

The explicit formula of the special function ^ , in the discrete field, is: 

/ IN 2fc + l 

o»l*- n <«] =( 1 + 3^ow) (19> 

where is always valid the following boundaru^l: 

j[k,n (i e -i)] -< 2 -< j[k,n(i e )] with o {h-i) = n (i t ) -i-,vi t , k,£e~N 



26 The letter (shin) is the thirteenth letter of the Arabian alphabet. 

27 The boundary can include the sign "=" if the integer variable k goes towards zero or the infinity. 
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The auxiliary integer function ft(7f), that really represents a growing "step function", 
is defined, for the intervals of 8 or 9 following values of k, in the following way: 

• ft (Ji) = for fc=l,...,8 

• ft (7 2 ) = 1 for fc=9,. . . ,16 ; ft (/ 3 )=2 for fc=17,. . . ,25 ; (7 4 )=3 for k=26,. . . ,34 

• n (7 5 ) = 4 for fc=35,. . . ,43 ; ft (7 6 ) =5 for k=U,. . . ,51 ; ft (7 7 )=6 for k=52,. . . ,60 

• ft (I 8 ) = 7 for fc=61,. . . ,69; ft (7 9 ) =8 for k=70,. . . ,78 ; O (7i )=9 for k=79,. . . ,86 

• ft (I u ) = 10 for k=87,. . . ,95 ; ft (Ii 2 )=H for fc=96,. . . ,104. ; etc. 

Successively we give (Fig. 1) the graphs, related to the families of ji functions, that are 
ft (7^)] and ji[A;, ft (7^ — 1)] , or better, to the set of the arcs belonging to them, and 
to the auxiliary function ft (7g) . 
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Fig. 1. 

Let's extend the dependence of the integer step function ft (7;) to the real field and let's 
use the following definition: 



ft (x) = min{/c G N: S fe+ i (x) > 2 } ; x G R + and where S k (x) = ( 1 + 

\ OX K ~\~ 1 



2x+l 
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By simple algebraic passages we have, by a more appropriate notation, due to Iverson, 



that: 



3x — 



1 



22^+T _ 1 

where \ x] means the smallest integer, greater than x or equal to it. 



(20) 



Therefore the special function jS possesses the following explicit formula in the real 



field: 



(x) 



1 + 



3x — Q (x) 



2x+l 



1 + 



V 



3x 



3x 



l 

i — 



\ 2x+l 
) 



Extending the field of definition of the variable k to the real positive numbers, it's 
possible to notice that such function, being represented by the union of continuous arcs 
(all above the straight line of height 2, see Fig.l) is actually assimilable to a piecewise 
continuous function. 

That being stated, an important subclass of completely monotonic functions consists of 
the Stieltjes transforms defined as the class of functions /: (0, oo) — > R of the form: 



1 1 \ [ d V (*) 

=a + J—t 



(21) 



where a > and fi (t) is a nonnegative measure on [0, oo ) with / < oo, see [2]. 



l+t 



In the Addenda and Problems in ([1], p. 127), it is stated that if a function / is holo- 
morphic in the cut plane A = C\ (— oo,0] and satisfies the following conditions : 

(i) 9/ (z) < for 9 (z) >- 

(ii) f(x) > for x >- 
then / is a Stieltjes transform. 
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2. THE REPRESENTATION AS A STIELTJES TRANSFORM 



In ([8], §6) the author characterizes the holomorphy (piecewise analytic) of the special 
function ^(z) in the cut plane B = C\ —73, —73 and proves a remarkable result that 
implies : 



lim J,{z) = 2 (zGB) 

x|^oo 



(22) 



To prove that the harmonic function satisfies $5^(21) < for 9 (z) >~ , we 

use that maximum principle for subharmonic functions, that can be found in ([4], p. 20), 
and show that lim sup of at all boundary points including infinity is less than or 

equal to 0. 



From (|22p we conclude that this is true at infinity. 

How, for definition ^(x) >~ for x >~ 0; these last statements imply the result (|2~T1) . 
The constant a in (|21 1) is given by : 



a = lim , i (2;) 

and therefore for the fundamental theorem of the special function ^2 we have, see ([8], §2): 

a = lim . Jb(x) = 2 
In (f2~T1) /i (i) is the limit in the vague topology of measures 

dfi(t)= lim -~$Sf(-t + iy)dt 

y^0+ IT 

For z £ B = C\ —73, — 73 we have in the close interval [—1,0] : 



rO) 



1 \ 2z+1 / / 1 

1 + ^— r ] = exp (2z + 1) • Log 1 + 



V ' 32 + 1 

where Log denotes the principal branch of the logarithm. 



3z + 1 



(23) 
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Let t G R and z £ C with 9 s (z) y 0. 

If z tends to t , then for ((TSJ) and (JS3J), results ( with £ G Z ^ ' 2lS • 

2t+l 



0" 



1 + 



cxp 



2 if t = 

(2t + 1) • log 



iftyO 



3f+2 
3i+l 



- k-in(2t + 1) toita fc = 1 i/ 
and with k = i/ — 1 ^ t ^ — | or i/ — | -< t -< 
2 if t = -l 



1 + 



3t-n 



2t+l 



In particular then we obtain , if y tends to + , foCj t G R: 
f if t < I/3 or t > 2 / 3 



1 



7T 



3/(-t + iy) 



! ((3f-l) 2 )' . 



(27Ti) • { 



2?r ((3<-2) 2 ) 

and using the identity (the Euler reflection formula): 

r(a) -r(l-a) 



3t-2 



3t-l 



3t^2 
3t- 



1} 



7T 



sin (era - ) 



we are now in a position to determine the following nearly exact integral representation ( 
Stieltjes transform): 



1 



3t - 1) 



2\ t 



r(2-t)-r(i-2-t) 



((3t 



St -2 



3t - 1 



3t-2 



3t- 1 



eft 



(x + t) 



(24) 



The approximation is essentially origined by neglecting the point of discontinuities of 
the first kind of the special function ^ , in the real field, between an interval and the 
following as far as the interval / [0, oo ) . 



28 Z denotes the relative integer set. 

29 In the discontinuity points 1/3 and 2/3 we respectively compute the limits of the real variable 



t on the left and on the right (see Fig. 2). 
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Successively we give (Fig. 2: x — > t) the graphs (red color) related to the function 30 ! 
M (t) in the real interval I [—1, 1]: for t = ^-fe I 1 (t) = and this point is a flex point with 
oblique tangent. 
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+ 
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V' 1 
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Fig. 2: the graph of ^ (t) 



That being stated, we denote the set of completely monotonic functions with C. 

Now, we also recall that a function /: ]0,oo[ — ► ]0,oo[ is said to be logarithmically 
completely monotonic [5] , if it is C°° and 

(-l) fc -[log/(^)] (fc) >0 for k = l, 2,3, ... 

To simplify we denote the class of logarithmically completely monotonic functions by C 
and the set of Stieltjes transforms by S . 

In order to prove that the special function ^J>(x) is logarithmically completely mono- 
tonic, we need the following lemma : 

S \ {0} C C 



30 £ ^ = M*l 
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This lemma is a consequence of the following result, established by Horn [6], that allows 
also to characterize the class of logarithmically completely functions as the infinitely divisible 
completely monotonic functions: 

Theorem 1: For a function f: ]0, oof — > ]0, oo[ the following are equivalent: 

(i) / G £ ; (ii) f a G C for all a y and a G R ; (iii) y 7 / G C for all n = 1, 2, 3... 

In fact, let/E<S(<ScC) and non-zero and let a >- 0, by Theorem 1 it is immediate 
to prove that f a G C. 

Now, writing a = n + a with n = 0, 1, 2, ... and < a -< 1 we have f a = f n ■ f a , and 
using the stability of C under multiplication and that / a G S =>• 5 \ {0} C £. 

In conclusion for (I24j) also the special function <ji(x) G £. 



3. THE CLASS OF BERNSTEIN FUNCTIONS 

There is an important relation between the set S of Stieltjes transforms and the class 
B of Bernstein functions. 

We recall that a function / : (0,oo) — > [0, oo) is called a Bernstein function, if / has 
derivatives of all orders and /' is completely monotonic. 

Now, if / is non-zero Stieltjes transform, then 1/f is a Bernstein function ([3], Prop. 
1.3). 

The special function G S \ {0} and this fact implies that is a Bernstein 

function. 

In addition, using the identity : 

1 /j ( x ) = X / x ■ u* ( x ) 

and remembering the following definition: 

A Bernstein function <j) is called a special Bernstein function if the function ^/^ ry\ is 
also a Bernstein function. 

we can conclude that x ■ ^ (x) is a special Bernstein function. 

The family of special Bernstein functions is very large, and it contains in particular the 
family of complete Bernstein functions (also known as operator-monotone functions, see [7], 
for instance). 
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Recall that a function (j) : (0,oo) — ^ Jft is called a complete Bernstein function if there 
exits a Bernstein function ij such that : 

0(A) = A 2 L [r? (A)] , A^O 
where L stands for the Laplace transform. 

Now, using the main results about the special function ^(x) ([8], §5 and §6) it is 
immediate to establish that x ■ ^J, (x) is a complete Bernstein function. 

Note also that a function / (x) is called a complete Bernstein function if, and only if, 



/ (x) = a + bx + J 



x 



t + x 



p(dt) 



(25) 



where a, b > and p is a Radon measure on (0, 00) such that J (1/(1 + i))p ((it) < 00 . 

0+ 

From this one, we may deduce that the function x — > / (x)/x is a Stieltjes transform 

[2]. 

This result was actually already obtained with the representation (flHj) of Stieltjes of the 
special function ji. 

At the end, recall that the following conditions are equivalent: 

(i) <p is a complete Bernstein function; 

(ii) (_\) is a complete Bernstein function. 

This result implies also that the function l/ ( ji(x) is a complete Bernstein function and, 



and l/jif^J are 



remembering the standard form (|25p and that the functions 1/ x ■ (x) 

Stieltjes transforms, it is easy and immediate to estabilish that the constant a (killing rate) 
and b (drift coefficient) are given by: 

= J!^Vu*f~) =Vu*(o) =~ 



6 = lim 1/ 







and the following representation formula which is analogous to the Levy-Khinchin formula: 

(3t - 2) 2 V r 3 t - 1 



i/us(z) 



1 1 

2 + 2 7 ] T(2-i) -r(l-2-t) 





((st-ify 



3t-2 



3t - 1 
3t-2 



x ■ t 

(t + x) 



dt 



34 



For the interplay between complete Bernstein functions and Stieltjes transforms we refer 
also to [9]. 

Finally, with an Euler-Venn diagram, we give the most important analytic properties of 
the special function ji . 




Completely Monotonic Functions ~ vs ~ Bernstein Functions 
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Abstract. We prove the following exact symbolic formula of the special function ^5, in 

the entire s-complex plane with the negative real axis (including the origin) removed, with 
a double Laplace transform: 

j(s) = L{2.6(t)+L{±-[j(t-e-™)- c i(t-e™)}}} 
where 5 (t) stands for the distribution of Dirac and e represents the Euler's number. 

1. THE EXACT SYMBOLIC FORMULA 

In [4] and [5] the author introduces and characterizes the discrete and special function 
j2 in the complex field. 

We recall that the explicit formula of the special function ^ , in the discrete field, is: 

1 \ 2k+l 



where it is always valid the following boundaru^\: 

j[k,n (i t - 1)] -< 2 -< j[k,n(i e )] with o {h-i) = n (i e ) -i-,vi t , men 
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The boundary can include the sign "=" if the integer variable k goes towards zero or the infinity. 
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The auxiliary integer function £l(Ie), that really represents a growing "step function", 
is defined, for the intervals of 8 or 9 following values of k, in the following way: 

• n (h) = for fe=l,...,8 

• n (I 2 ) = 1 for k=9,. . . ,16 ; Q (J 3 )=2 for k=17,. . . ,25 ; ft (J 4 )=3 for k=26,. . . ,34 

• n (J B ) = 4 for fe=35,. . . ,43 ; ft (I fl ) =5 for fc=44,. . . ,51 ; (J 7 )=6 for fe=52,. . . ,60 

• (J 8 ) = 7 for jfe=61,. . . ,69; ft (J 9 ) =8 for fc=70,. . . ,78 ; 17 (Ji )=9 for jfe=79,. . . ,86 

• ft (In) = 10 for k=87,. . . ,95 ; O (Ji 2 )=ll for jfe=96,. . . ,104. ; etc. 



In particular the author proves in ([4], §§5-6) that the special function ^J, in the real 
field is completely monotonic for all the closed intervals It and therefore the following first 
approximation "~" in the complex field: 



J (s) « L s [F (t)} =L S [2-Q (t)] = J e~ st d [2 ■ 9 (t)} = 2 



(26) 



where L s denotes the Laplace-Stieltjes transform and [t] represents the unit step function 
or Heaviside's function ( 9 [ t < ] = ; 8 [ t > ] = 1 ) 

Now, we observe that the Laplace-Stieltjes transform is closely to related other integral 
transforms, including the Fourier transform and the Laplace transform. 

In particulary if g has derivative g', then the Laplace-Stieltjes transform of g is the 
Laplace transform of g'. 

Consequently, considering the derivative of Heaviside step function, from (|26p we have : 



ji (a) w L [2 ■ 8 (t)] = J e~ st 25 (t) dt = 2 



(27) 



where L denotes the ordinary Laplace transform and 5 (t) stands for distribution of Dirac. 
In ([5], §2) the author find practically the following nearly exact integral representation: 



3t - 1)' 



2^ t r 



r(2-t).r(i-2.t) ( f3t _ 2 



3t-2 



3t- 1 



3t- 2 
3* - 1 



(x + t) 



(28) 



where T denotes the Eulerian gamma function (or the Eulerian integral of second kind). 
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The approximation is essentially origined by neglecting the point of discontinuities of 
the first kind of the special function ^ , in the real field, between an interval 1% and the 
following as far as the interval / [0, oo ) . 

We recall that extending the field of definition of the integer variable k of the discrete 
and special function ji to the real positive number, it's possible to notice that such function 
is actually assimilable to a piecewise continuous function. 

On the contrary of equation (|27p . in the equation (|28|) there is an Stieltjes transform, 
that arises naturally as an iteration of the ordinary Laplace transform. 

In fact, if 

f( x )= yv*v (t)dt 



where 

oo 

ip (x) = / e~ xt i) (t) dt 



o 



then, changing the order of integrations in the double integral by appealing to Fubini's 
Theorem, we have formally ([1], p. 127 and [8], p. 335): 

oo oo 

f(x)=L{L [if>t]} = J e~ xu du J e~ ut ^(t) dt (29) 



i/i (t) dt J e ~ u{x+t) du 





Hence 



f(x)=[^dt (30) 
J x + t 


This last equation we refer to as Stieltjes transform, or from another point view as the 
Stieltjes integral equation. 

However, we shall usually be concerned with the more general case in which the integral 
equation ([30]) is replaced by a Stieltjes integral: 







da (t) 
x + 1 



and in this form the equation was considered by T.J. Stieltjes in connection with his work 
on continued fractions [6]. 
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In (|28l) the Stieltjes integral, that is analogous to (|30p . converges and then we have (see 
[8], Theorem 7b, p. 340): 

lim f(-t-iri)-f(-Z + iri) = ^jj+V^) m 

for any positive £ at which tp (£+) and i/> (£— ) exist. 
For, simple computation gives : 

oo 

f(-t-iv)-f(~t + iv) 1 / 



2vr i tt J (t - £Y + r/ 2 



(it 



o 



The integral is known as Poisson's integral for the half-plane or as Cauchy's singular 
integral ([7], p. 30). 



1> (x) (32) 



The result (|3ip may also be written symbolically as 

/ (x ■ e~ i7r ) - f (x ■ e m ) 

In fact it is sufficient to compare the equation 11.8.4 in ([7], p. 318). 

That being stated in ([4], §6) the author characterizes the holomorphy (piecewise 
analytic) of the special function j*(s) in the cut plane A = C\ —73, — 73 • 

In the close interval [—1,0] results: 

2s+l 



<»» = (' + ^Tl) ' ^ " exp ( (2s + 1} ' L " a { l + 3JTl)) (33> 



where Log denotes the principal branch of the logarithm in the interest of obtaining a single 
value mapping. 

In effects the special function possesses two branch points at s = — | and s = — | 
and therefore a cut, joining the two branch points, would prevent s to circulate around 
them, and the special function j3 (s) can be treated as a holomorphic function. 

We observe that, in our particular case, the results (|3"2"j) and (|3"T]) with (f3"3"j) can be 
written and compute in the following way: 

Ji (x ■ e~ ln ) - JS (x ■ e m ) Ji (-z - iy) - j5 (-a + iy) 

= lim 

2m y^o+ 2m 
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_ r f exp((2(-x-»y)+l)-Log(l+ g )+1 )) -exp((2(-x+iy)+l)-L 0g (l+ 3( _^ iy)+1 )) 

- lim^ + | 2^ 

sin (2vrf) (g* ~ if) ' f 
2 ^ ((3i - 2) 2 ) * I 

Finally, with ([2T ]) . (pg|hP^ .([M l) and using the identity (the Euler reflection formula): 

r (a) • r (1 - a) = 

sm (air) 

we can prove, for the linearity property of the Laplace transformation, the following exact 
symbolic formula of the special function j*, in the cut plane A = C\ (— oo,0]: 

j(s) = L{2.6(t) + L{±- i .[j(t-e-™)-j(t-e™)]}} (35) 

where S (t) stands for the distribution of Dirac and e represents the Euler's number. 

In conclusion we give also the three-dimensional graph of absolute value of the special 
function ^J, of a complex variable s = x + iy (Fig. 1) from DERIVE computer algebra 

system and in APPENDIX we explain the choice of the Arabian letter ^ . 



3t-2 
3t-l 



3t-2 
3t - 1 



(34) 




Fig. 1 
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2. ADDITIONAL ANALYTIC REMARKS 



The beautiful symbolic formula (35) is in practice a consequence of the fact that Laplace- 
Stieltjes transform are often written as ordinary Laplace transform involving the distribution 
of Dirac, sometimes referred to as the symbolic impulse function 5 (t). 

Besides that the construction-algorithm ([4], §1) of the discrete and special function ^ 
had already evidenced the necessity of an asymmetric impulsive force in order to maintain 
the set of the points belonging to them all above the horizontal straight line of height 2. 

Now, we observe that the asymmetrical impulse function (5+ (i) is more suitable for 
use in connection with the one-sided Laplace transformation that the symmetrical impulse 
function 5 (t). 

In effects in the formula (35) the presence of the symbolic impulse function (5+ (t) would 
be more appropriated , also because if one applies the Laplace transformation to the "defi- 
nition" of the impulse function 5 + (i), one obtains the formal result: 



We recall that the asymmetrical impulse function 5+ (t) is defined ([3], see Sec. 21.9-6) 



L[S+(t)] = 1 



by: 



{ 





/(x+0) 



if x^ia or x>b 
if a<x-<b 



{a^b) 



It is possible to write: 



5+ (t) 




where U+ (x) denotes the asymmetrical unit-step function : 




Successively we give the cartesian diagram (Fig. 2) of the unit-step function [/+ (x) and 
of the one related to approximation of the impulse function <5 + (t). 
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In the papers [4], [5] and in this paper it has emerged all the importance of the tools 
like the Laplace transformation and the Stieltjes integral with certain facts from the theory 
of functions of a complex variable in order to characterizing the analytic properties of the 
completely monotonic functions. 

In particular we recall that the Stieltjes integral is much used in Mechanics and Proba- 
bility, since it unifies the treatment of the continuous and discrete (and mixed) distributions 
of mass or probability. 

If a (x) is piecewise differentiable, then da (x) = a' (x) dx, and the Stieltjes integral is 
simply in the following form (reduction of a Stieltjes integral to a Riemann integral): 

b 

J f (x) a' (x) dx 

a 

For instance a real and similar case to ours is the following: if a (x) is a Heaviside step 
function, with point masses rrn at x = Xi, then 

da {xi) = lim [a (xj + e) — a (xj — e)] = rrii, 
In this case the integration by parts is usual: 




b b 
J f{x)dg{x) = f(x)-g(x)\ b a J g (x) df (x) 

a a 
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Suppose a (0) = , a(x) = o (e cx ) as x — > oo , and that Sis > c, then 

oo oo 

J e~ sx da (x) = s- J a (x) e~ sx dx (36) 



The integral on the left side represents a Laplace-Stieltjes transform, while the integral 
on the right side is an ordinary Laplace transform. 

More precisely (see [4], compare the equation 14) the second member of (I36p is called 
s-multiplied Laplace transform. 

However, we also recall the following result, due to S. Bernstein ([2], pp. 439-440), who 
was the starting point of many researches in the Probability theory: 

Theorem: A function ip (s) on (0, oo) is the Laplace transform of a probability distri- 
bution F (x): 

oo 

il){s) = j e~ sx dF(x) 
o 

if and only if it is completely monotone in (0, oo) with ip (0+) = 1. 

In particular, an immediate example for this beautiful theorem is the following function: 




Notice that the Dirac delta function may be interpreted as a probability density function 
and that the cumulative distribution function is the Heaviside step function. 

It is known that if X is a random variable, the corresponding probability distribution 
assigns to the interval [a, b] the probability Pr[a < X < b]; for example the probability that 
the variable X will take a value in the interval [a, b]. 

Now, the probability distribution of the variable X can be uniquely described by its 
cumulative distribution function F(x), which is defined by: 

F (x) = Pr [X < x] 

for any ieR and where the right-hand side represent the probability that the variable X 
takes on a value less than or equal to x. 



We observe that, in our case, the Heaviside step function is the cumulative distribution 
function of a random variable which is almost surely 0. 



44 



APPENDIX 

The choice of the Arabian letter is simply a consequence of the following four reasons: 

• the exhaustion of Latin, Greek, Gothic, Jewish and other letters to name a new special 
function; 

• the reference to the strokes of continuous curves (small arcs) that characterize the 
same function ( [4], see Fig. 3); 

• the three dots over the letter that one by one remember the three second-order Eu- 
lerian numbers: 1, 120, 494, inside the Eulerian triangle ( [4], see Fig. 6); 

• the meaning of such a letter, that in Al Karaji's algebra (about 1000 A.D.j^l was 
declared as the unknown "par excellence" , that is absolutely comparable to our "x" , 
but also with the characteristic to form the only conjunction- "ring" between the world 
of the unknown (algebra) and the world of the known (arithmetics). 

Actually it seems that the origin of the symbol ^ , that is pronounced shin , is in the 
Ancient Egypt. 

Mi! 

In fact the hieroglyph td=! = t = l 1 that is the Egyptian syllable sha, is similar and it represents 
papyrus plants along the Nile. 

Acknowledgements. This paper and the papers [4], [5] are dedicated to the memory 
of my father Alessandro Ossicini. 
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